Instability and spatiotemporal rheochaos in a shear-thickening fluid model 
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We model a shear-thickening fluid that combines a tendency to form inhomogeneous, shear-banded 
flows with a slow relaxational dynamics for fluid microstructure. The interplay between these factors 
gives rich dynamics, with periodic regimes (oscillating bands, travelling bands, and more complex 
oscillations) and spatiotemporal rheochaos. These phenomena, arising from constitutive nonlinearity 
not inertia, can occur even when the steady-state flow curve is monotonic. Our model also shows 
rheochaos in a low- dimensional truncation where sharply defined shear bands cannot form. 



Complex fluids exhibit much interesting behavior un- 
der shear, due to strong couplings between mesoscopic 
structure and flow. Experiments show cases where, un- 
der steady external driving, an unstable flow arises, giv- 
ing a time-dependent strain rate at constant imposed 
shear stress, or vice versa. Sustained temporal os- 
cillations are seen in surfactant mesophases and solu- 
tions 0, 0, El & El and polymer solutions 7], while 
erratic temporal responses have been found both in these 
and in related materi als, e.g., wormlike micelles pllfl Ho|. 
lamellar phases llMEl and colloids [13,0. There are 
strong indications that these erratic signals result 

from a deterministic chaotic dynamics. Chaotic behav- 
ior of bulk flows at virtually zero Reynolds number (neg- 
ligible inertia) must stem from nonlinearity within the 
rheological constitutive equation, and has been dubbed 
'rheochaos' BE! El- 

Such flow instabilities affect both shear-thinning || 
and shear-thickening micellar materials Shear-thick- 
ening is, in itself, a widely observed but poorly under- 
stood phenomenon which affects not only micelles but, 
e.g., dense colloids, where there is again evidence of bulk 
rheological instability 0, El • Below we study a simple 
model for a shear-thickening fluid at steady controlled 
stress, which generalises to spatially inhomogeneous flows 
a model first proposed in Ref. [l4j . The latter was shown 
to give simple oscillations but not chaos (unless an uncon- 
vincing 'double memory' term was used). We find that 
the interplay of a very simple structural memory with 
constitutive nonlinearity can, as hoped by the authors 
of [3 , give complex dynamics and rheochaos - but only 
if spatial heterogeneity is allowed for. Our work is related 
to, but different from, that of Fielding and Olmsted [H| 
which addresses shear-thinning fluids 17]. Also related 
is Ref. [13 which concerns nematic liquid crystals (which 
can be chaotic even without spatial heterogeneity). We 
show below that there are robust generic features in the 
rheochaos produced in these various models, which gen- 
erally involve a failed attempt to create a steady shear- 
banded flow. (In the shear thinning case the bands have 
a common stress but unequal strain rates; for us, the 
reverse is true.) One distinctive finding of our work is 
rheochaos in a model where the steady state flow curve is 
strictly monotonic. Here (in contrast with [l£|) no steady 



banded solution, stable or otherwise, exists. Yet an in- 
nate tendency to form bands transiently, combined with 
structural memory that frustrates their long-term sur- 
vival, can lead to chaos. A second unexpected finding is 
that a low dimensional truncation of our model also gives 
chaos, even though sharp interfaces between shear bands 
are then suppressed. Thus high-order Fourier compo- 
nents of spatial variation are not essential to rheochaos. 

Candidates for a slow structural mode include the 
mean length of wormlike micelles, local composition vari- 
ables (e.g., colloidal volume fraction) and 'fluidity' pa- 
rameters [l^l LLil reflecting local microstructure, bonding 
state, etc.. In all cases, the time scale rs is distinct from 
the Maxwell time tm for linear stress relaxation; and in 
our work, we assume ts/tm > 1 (in contrast to 16]). An 
involvement of slowly evolving fluid structure has been 
clearly evidenced in some experimental cases 0. CI Q . 

The model. We assume that the shear stress a de- 
couples from other stress components 20] and depends 
only on the rate of shear strain 7. For a cylindrical 
Couette geometry with axial coordinate z, we consider 
one-dimensional inhomogeneity along this direction only 
(as would arise for steady shear bands in a thickening 
material 21]). We assume the shear flow is homoge- 
neous within each slice of height z, appropriate to a low- 
Reynolds limit; 7 is then uniform and fixed by the im- 
posed wall velocity. The shear stress cr(z,t) evolves as 
follows, where t is time and units are such that the tran- 
sient elastic modulus is unity: 

&(z,t) = j-R(a)-X [ M(t-t')o-(z,t')dt' + KV 2 o- (1) 



The term R(cr) = aa — bcr 2 + ccr 3 corresponds to instan- 
taneous nonlinear stress relaxation [l4j, with a = 1/tm, 
and 6, c chosen so that cr(R) is an S-shape (see Fig. 1). 
This creates a tendency for the fluid to form shear bands 
stacked in the z-direction. The integral term represents 
a retarded relaxation of stress and/or a slow evolution of 
fluid structure on timescales rs (see [lj)- .M is a decay- 
ing memory kernel, chosen as Ai(t) = r^ 1 exp(— t/rs); 
A > governs its strength. Finally, the nonlocal term, 
assigning to stress a diffusivity ft, selects a unique banded 
flow in steady state [53| . Note that our model is closely 
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FIG. 1: Long-term flow curve R(a) = 7 — Act (thick line), and 
underlying short-term curves obeying R(a) — 7 — Am (thin 
lines, with, from left to right, ra = 0,l,2,..., 20). The stress 
range between the dotted lines corresponds to the unstable 
region (for rs = 100). Parameters: A = 40, a = 1/tm = 100, 
6 = 20, c= 1.02, « = 0.01, H = 1. 



analogous to the FitzHugh-Nagumo for neuronal activ- 
ity [23|, albeit with an unusual long-range coupling. 

Qualitative features. We now rewrite equation JJJ as 
an exactly equivalent differential system: 
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FIG. 2: Phase diagram of the model when rs and (a) are var- 
ied: (A) chaotic point, (•) periodic point. Three main regimes 
are observed: (O) oscillating shear bands, (T) travelling shear 
bands, (C) chaotic regions. The outer dashed line is the linear 
stability limit, R'((a)) + l/r s = Ql|2f]. The dotted lines 
between regions are guides to the eye; there are no sharp tran- 
sitions. The C-regions enclose all observed chaotic points, but 
contain internal structure with periodic and chaotic pockets. 
Numerical parameters as in Fig. 1. 



& = 7 — R((j) — Am + ftV 2 a and rh = (2) 

rs 

The integral term in eq. has become an auxiliary vari- 
able m(z,t), which we call the memory. The memory 
encodes the delayed (structural) part of the stress relax- 
ation, and, for our choice of decays with rate r^T 1 to- 
wards the local value of the stress a(z, t) Q. The homo- 
geneous version of the present model (no ^-dependence) 
shows temporal instability, leading to van der Pol type 
oscillations Q. 

The steady-state flow curve a(j) for the model is ob- 
tained from a = r'n = V 2 a = 0, as R(cr) + Xa = 7. At 
each point on the curve, the memory m has relaxed to 
the steady state stress (m = a). However, this relaxation 
involves the long timescale rs. At times much shorter 
than rs , the fluid will instead behave as if the memory m 
were frozen: one has a set of 'instantaneous' flow curves 
R((j) + Am = 7. Thus, despite the fact the steady-state 
flow curve is monotonic, the existence of fixed m curves 
with an S-shape gives rise, on short times, to spatial in- 
homogeneity in the form of shear bands along z (Fig. 1). 
As these bands are subjected to the local van der Pol 
instability referred to above, a very rich spatio-temporal 
dynamics arises as we now show. 

Fourier- Galerkin truncation. Equations @ were 
solved numerically using a spectral Galerkin trun- 
cation |23, with a(z,i) and m(z,t) decomposed 
in Fourier modes cr^.(t) and rrik(t): o~(z,t) = 
^2k=o a k(t) cos(k7rz/H) etc., with H the axial extent of 
our Couette device (0 < z < H) and N the order of 



the truncation. (Sine modes are excluded by zero stress- 
flux boundary conditions, V<r = 0.) The truncation thus 
replaces eqs. for a(z, t) and 771(2, i) by 2N coupled or- 
dinary differential equations for the modes cr^ (t) , rrik (i) . 
The results presented below used a high-resolution trun- 
cation (N = 40); smaller 7V(> 20) already give simi- 
lar results. We work at fixed torque on the Couette, 
that is, at fixed value of the spatially averaged stress 
(a)=J H a(z,t)dz/H 0- 

Phase diagram. In Fig. 2 we show a 'phase diagram' 
grouping the different types of spatiotemporal dynamics 
seen 27] on varying (a) and ts/tm (with tm = 0.01, and 
remaining parameters held at the values of Fig. 1, and 
small amplitude noise in the initial condition). Remark- 
ably, despite the high-dimensional dynamical system un- 
der consideration (2N = 80), the phase diagram displays 
a simple overall structure, with three main regimes: pe- 
riodic response with oscillating shear bands at extremely 
long rs; periodic response with travelling bands at long 
rs; and chaotic response at shorter rs and off-centered 
values of (a). 

Oscillating shear bands. The succession of behaviors 
observed along the horizontal line ts/tm = 10 4 in the 
phase diagram typifies this regime. Near the middle of 
the line (e.g. (a) — 7.0) we find two 'flip-flopping' bands: 
as Fig 3-a shows, mid-cycle a low-shear band and a high- 
shear band each span about half of the cell. Each under- 
goes a local van der Pol type oscillation; to keep (a) fixed 
as prescribed, these are synchronous and each half-cycle 
the identity of the two bands is reversed. The time series 
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FIG. 3: Typical responses in the oscillating shear band 
regime, chosen along the line ts/tm = 10 4 in the phase dia- 
gram, (a) Flip-flopping bands at (a) — 7; (b) Zig-zagging in- 
terface at (a) =7.1; (c) Complex periodic motion at (a) = 9. 
Each group presents time series of the stress a at z = 2/3, 
the shear rate 7, and a space-time plot of a(z,t) with t verti- 
cal, z horizontal (clear shades correspond to high stress, dark 
shades to low stress). 



of the local stress (e.g., <r(2/3,£); Fig 3-a) shows a large 
amplitude oscillation, close to a square wave, with a flip- 
flop period rp of order rs. The shear rate 7 (independent 
of z) accordingly shows a periodic evolution, with a more 
complicated waveform. 

If we now raise or lower (a) slightly (e.g., (a) =7.1, 
Fig. 3-b), in addition to the bands' flip-flopping, their in- 
terface adopts a zig-zag motion. (Imposition of (a) now 
enforces unbalanced proportions of the low- and high- 
shear bands; hence on flipping, the interface must move 
to and fro to maintain these proportions throughout the 
cycle.) The time series <r(2/3, t) is slightly distorted, with 
period now a multiple (here three) of rp. Moving further 
towards the wings of the phase diagram, the intrinsic 
flipping dynamics of the bands, coupled to the global 
constraint on (a) (e.g., (a) = 9.0, Fig. 3-c) gives an ex- 



FIG. 4: Time series and space-time plot of the stress a(z,t) in 
the travelling band regime. Parameters: ts/tm = 90, (a) = 7, 
others as in Fig. 1. 
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: Time series and space-time plot of the stress cr(z,t) 
chaotic bands regime ('bubbly' type). Parameters: 
= 90, (a) = 3.6, others as in Fig. 1. 



tremely complex motion which, remarkably, manages to 
remain periodic. The time series <j(2/3,t) becomes very 
ragged and the period is a large multiple of rp. 

Travelling shear bands. For shorter ts (20 < ts/tm < 
10 3 ) we find a new periodic regime in which shear bands 
nucleate at a boundary and then cross the system, usu- 
ally one at a time, at roughly constant velocity (Fig. 4). 
Rarely, bands can nucleate periodically from a non- 
boundary point 26], giving two outgoing bands in a 
one-dimensional analogue of the 'target patterns' seen 
in chemical oscillators |28| . Details here depend strongly 
on the initial noise. The travelling bands are kine- 
matic waves arising from staggered phase distribution 
of the local van der Pol-like oscillation. Our data [2(| 
also indicate that the band velocity decreases as ts is in- 
creased. Somewhat similar band motion was reported in 
Ref.|lq but there the bands 'ricochet' off the walls of the 
container; ours disappear and renucleate at the opposite 
wall. 

Spatiotemporal rheochaos. The third regime is chaotic 
and arises mainly for 10 < ts/tm < 100, in the wings of 
the phase diagram. A typical dynamics in this regime 
(Fig. 5) comprises a 'bubbly' phase of localized, short- 
lived shear bands that appear and disappear erratically 
in the cell. (Occasionally a band survives longer and 
shoots across the cell.) Other types of chaotic space-time 



4 



patterns were also found [2f|: 'wiggling' travelling bands 
and 'defect dominated' regimes (the latter close to those 
found for shear thinning micelles Ref. 16]). Chaos, for 
each point marked in Fig. 2, was confirmed by a positive 
Lyapunov exponent [2y |. 

Role of stress diffusion. The diffusivity n has a singu- 
lar role in steady-state shear banding, selecting a unique 
banded flow |22|. However in our work stress diffusion 
is, in Fourier space, a diagonal term acting mainly by 
damping the higher modes: small n increases the order 
N required for realistic truncation, and this tends to favor 
chaos. We suspect that n also plays a role in the crossover 
from travelling to oscillating bands seen at ts/tm — 10 3 
in Fig. 2: for our parameters, this happens roughly when 
rs becomes equal to the typical diffusion time across the 
Couette, Tdiff = 1/nq 2 (with q = tt/H). We do find the 
crossover to move towards longer timescales as k is re- 
duced [26j |. 

Low-truncation limit; route to chaos. Our work, with 
that of [01 E3 suggests a generic interpretation of 
rheochaos in terms of the erratic motion of discrete band 
interfaces But if sharp interfaces were essential, 

there would be no chaos in a low-mode truncation that al- 
lows only smooth variation of a, m on the scale of the sys- 
tem. To test this, we studied the case N = 3; with o~o and 
mo fixed by (<r), our dynamical variables comprise <Ti(t), 
o~2(t), Trti(t) and ra2 (£). Our numerics show 26] that, 
although shrunk, the chaotic regions persist. Sharply re- 
solved band interfaces are thus not essential to rheochaos. 
Unlike the (numerically more delicate) high-resolution 
truncation, the N = 3 case allows determination of the 
transition to chaos as (a) is varied at fixed rs: we found 
a classical period-doubling cascade. 

Concluding remarks. Our minimal model of a shear- 
thickening fluid combines the tendency to form shear 
bands with structural relaxation. This combination leads 
to oscillatory and travelling banded solutions and to spa- 
tiotemporal chaos. Our work strongly supports a generic 
'frustrated shear-banding' picture of rheochaos (e.g. LLJ) 
as instanced in recent work on two other models |l5lllq: 
it seems increasingly clear that spatial inhomogeneity 
plays a key role in rheochaos. Intriguingly, in our model 
it is enough to represent this by a very low dimensional 
truncation of the spatial structure: transient shear bands, 
in the sense of sharp interfaces between layers of mate- 
rial of different flow history, are not essential. But, in 
practice, the stress diffusivity n usually is small and en- 
sures rather sharp interfaces; we suspect the 'frustrated 
banding' picture to hold in many of the experimental re- 
alizations of rheochaos so far reported. 
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